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Abstract 

Based on the approach of flow distances, the international trade flow system is studied from 
the perspective of multi-layer flow network. A model of multi-layer flow network is proposed 
for modelling and analyzing multiple types of flows in flow systems. Then, flow distances are 
introduced, and symmetric minimum flow distance is presented. Subsequently, we discuss the 
establishment of the multi-layer flow networks of international trade from two coupled viewpoints, 
i.e., the viewpoint of commodity flow and that of money flow. Thus, the multi-layer flow networks 
of international trade is explored. First, trading “trophic levels” are adopted to depict positions 
that economies occupied in the flow network. We find that the distributions of trading “trophic 
levels” have the similar clustering pattern for different types of commodity, and there are some 
regularities between money flow network and commodity flow network. Second, we find that active 
and competitive countries trade a wide spectrum of products, while inactive and underdeveloped 
countries trade a limited variety of products. Besides, some abnormal countries import many types 
of goods, which the vast majority of countries do not need to import. It may indicate an abnormal 
economic status. Third, harmonic node centrality is proposed and we find the phenomenon of 
centrality stratification. It means that competitive countries tend to occupy the central positions 
in the trading of a large variety of commodities, while underdeveloped countries likely in the 
peripheral positions in the trading of their limited varieties of products. Fourth, we find that 
manufactured products have significant larger mean first-passage flow distances from the source to 
the sink than that of primary products. 
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I. INTRODUCTION 


Flow network is an important tool for describing and analyzing open flow systems and 
has been studied extensively in a large range of flow network systems, such as ecological flow 
networks MUD], world trade flow networks mi, traffic and transportation flow networks 
[T2J . Flow networks are commonly modelled by directed weighted networks, where directions 
and weights of edges represent directions and volume fluxes of flows respectively. Because 
open flow systems always exchange energy, matter and information with their surroundings, 
flow networks normally have two special nodes (i.e., the source and the sink) representing 
the environment, where all flows are supposed to be from the source, through the system 
and finally to the sink node [lllj. 

Based on the model of flow networks, many useful methods have been developed for 
exploring structures and dynamics of flow systems. For example, Rosvall and Bergstrom [7] 
proposed a method of probability flow of random walks for revealing community structure 
in weighted and directed networks. Vezina and Platt [8] described an inverse method for 
estimating network fluxes in undersampled environments. In a recent work [9], we solved 
a problem that how to measure the distances between nodes in flow systems. Several flow 
distances, i.e., the first-passage flow distance, the total flow distance and the symmetric 
flow distance, were put forward, and their helpfulness in calculating “trophic levels” and 
clustering for nodes has been preliminarily shown in the work. 

Using these methods, some significant laws and important knowledge have already been 
revealed in such kinds of complex flow systems. For instance, in living things, such as animals 
[lj, plants [2j] and microbes [3j, allometric scaling laws (e.g., the power-law relationship 
between an animal’s metabolic rate and its body mass P) exist extensively due to energy 
transportation on energy flow networks of living organisms M- The recent works have 
shown that the universal allometric scaling law also exists in a much broader range of flow 
networks, e.g., weighted food webs [9]J and trade flow networks HU 

Currently, most of studies on flow networks focus on monolayer flow networks 0EUMI]. 
However, complex systems containing multi-flows cannot be well represented by monolayer 
flow networks. For example, in a transportation system, aircraft traffic, train traffic and 
highway automobile traffic may coexist; in a city infrastructure, water flow, power flow and 
gas transmission flow exist at the same time; in a human body, the circulations of blood, 
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lymph and interstitial fluid coexist. If we do not distinguish between different types of flows, 
and simply adopt monolayer flow network to describe these flow network systems, many im¬ 
plicit information of multi-flows will be lost. Thus, many important knowledge for multi-flow 
systems may not be revealed. Besides flow systems having multi-flows of different types of 
substances, complex flow systems containing flows of the same substance with different labels 
(e.g., groups) also need a practical tool for system modelling and analytics. For instance, in 
a world wide trade system, money flows for trading different types of commodities may flow 
among different countries; in the logistics industry, different transportation flow networks 
are adopted by competing logistics enterprises. In such cases, multi-flows cannot be well 
distinguished with each other if we simply use monolayer flow network to model the systems. 
Overall, because these co-existing multi-flows are hard to be finely modelled by a monolayer 
flow network, it is necessary to introduce a new concept called Multi-layer Flow Networks 
(MFNs). Examples of MFNs include multi-layer trade flow networks, where different types 
of commodities flow in different layers, and multi-layer transportation flow networks, where 
passenger traffic flows on the layer of air transportation network and that of rail network. To 
the best of our knowledge, a generalized theoretical framework of multi-layer flow networks 
is still lack. So it is necessary to study multi-layer flow networks and their applications in 
depth. 

In this paper, we use world wide trade as a specific instance of multi-layer flow network. 
In previous work, international trade has been studied from the viewpoint of networking. 
Different types of trade networks (such as binary or weighted, directed or undirected ones) 
have been built to model real trade interactions between countries. For example, Serrano 
and Boguna [16] finds that the binary and undirected world trade web also presents the 
typical properties of complex networks, e.g., scale-free inhomogeneities and a high clustering 
coefficient. Fagiolo et al. S3 studied the topological properties and their evolution over time 
using a weighted world trade network. Fan et al. [T8] explored countries’ roles and positions 
using the improved bootstrap percolation and the other methods in an international trade 
network of single-layer. Unlike these previous studies HHHMH], here we use the framework 
of multi-layer flow network to explore novel characteristics of international trade networks 
based on the methods of flow distances, which is recently proposed in [l]3j. 

The rest of the paper is structured as follows. In Section 2, a formal description for 
multi-layer flow networks is presented, and flow distances are introduced. In Section 3, we 
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use multi-layer flow network of international trade as a case, and discuss the establishment 
of the multi-layer flow network from the viewpoints of trading commodity flow and money 
flow respectively. Then, using the approach of flow distances, some results on countries’ 
trophic levels, node centralities and mean first-passage flow distances from the source to the 
sink are discussed in Section 4. At last, we give the conclusions of this study in Section 5. 


II. METHODS 

A. Modeling multi-layer flow networks 

Multi-layer flow networks (MFN) can be regarded as a special type of multi-layer networks 
mm- This tool is able to well depict multiplex flows in different layers. Thus, we have 
the following definition. 

An multi-layer flow network is a pair MFN = (G, E ), where G = {G a \ a G 0,..., M — 1} 
is a family of directed (binary or weighted) graphs G a = { V, E a } (called layers of MFN), 
and 

E = {E af 3 C (v,v);v G V;a,f3 G {0, ••• ,M -l},a^ (3} (1) 

is the set of interconnections between nodes and their counterparts in the rest of layers. 
The node sets V in each layer are all the same, and are supposed to contain N common 
nodes and two special nodes “source” and “sink”, where “source” is the start of all flows 
and denoted as node 0, and “sink” is the end of all flows and labeled as node N + 1. E a , 
the adjacency matrix of layer a, can be written as follows. 

E a = {efj}( n+ 2 )x(n+ 2 )', b j G {0, • • • , N + 1}, a G {0, • • • , M — 1}. (2) 

where ei t j is the flow from node i to node j. Especially, the first column and the last row 
of this matrix are all 0 because there are no inflow to “source” node and no outflow from 
“sink” node. Normally eo,v+i is 0, because ordinarily there is no flow from “source” to 
“sink” directly. The total inflow to node j (denoted as e.j) is calculated as ^^Lq 1 and 
the total outflow from i (labeled as e*.) is ^2^=0 e *r Because the flow system is assumed to 
be in equilibrium, except “source” and “sink”, all other nodes have a balanced inflow and 
outflow, that is e. t = e*., where i — 1, • • ■ , N. The flows to “sink” (i.e. , e^v+i) are regarded 
as dissipation. 
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FIG. 1: Illustration of a multi-layer flow network 


B. Distance on multi-layer flow networks 

In [IS], L. Guo et al. have proposed flow distances on open flow networks. In this 
subsection, we will extend them to those on multi-layer flow networks. 

Let us consider a single layer cr of a given MFN firstly. Suppose a large enough number 
(say A) of particles flow along directed edges randomly. This random walk process can 
be assumed as a Markov chain. The mean first-passage flow distance (MFPFD) from one 
node i to another node j (denoted as kfl is defined as the expected number of steps for 
reaching j for the first time, given that initially the particles are at i. The mean total flow 
distance (MTFD) (denoted as tifl is the average number of steps for arriving at j regardless 
of whether it is the first time of arriving, also given that the particles are at i initially. 

To illustrate the above concepts vividly, the following thought experiment is designed 
and considered. Suppose all particles are initially no color. There are two situations. 1) If 
particles go through node i, their color turns blue; then when blue particles arrive at node 
j, they become red. In this situation, MFPFD from i to j is the average number of steps 
of blue particles turning red. 2) If particles go through node i, their color turns blue; blue 
particles remain blue, when arriving at node j. Then, MTFD from i to j is the average 
number of steps of all accumulated blue particles (also containing those whose arrival times 
is more than 1) arriving j. Here, a blue particle may be counted repeatedly if it arrives at 
j for multiple times. 
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Based on the above experimental description, the computation of flow distances are given 
as below. Formally, we have a stochastic matrix describing the transitions of the Markov 
chain M a = {m?-}( JV + 2 )(jv+ 2 ), where 


m ij = 


Y^ iV “I 

/ j j =1 


N+l CL 
3=0 C ij 


^N+l a 
/-^j=0 

S^N+1 a 

Z-/j =0 


7^0 
= 0 . 


(3) 


Then, for this absorbing Markov chain, its fundamental matrix is defined as below. 


U Q = (I- M a )~ l = I + M a + M 2 q + • • • 


(4) 


where / is the corresponding identity matrix with the same size of M a . Thus, we have the 
matrix of MTFD T a = {f"-}(jv +2 )(w+ 2 ) and that of MFPFD L a = {Z"}(jv+ 2 )(iv+ 2 ), where 

(mjjI y 

y U„ 

lOi _ J.CX. _ J.OL 

L ij — L ij L jj * 


(5) 

( 6 ) 


The detailed derivations of the above formulas can be found in [T3j. Because and /"• are 
asymmetric and cannot well satisfy the applications which need symmetric distance metrics 
(e.g., clustering and generating minimum spanning tree), a symmetric flow distance c“- was 
also introduced in [13], which is given as below. 


c 


a 

ij 


nia la 

a i3 l 3i 
/“• + /"• ’ 


(7) 


Here, we present an alternative symmetric flow distance called symmetric minimum flow 
distance (SMFD), which is simply calculated as below. 


fij = l %}' ( 8 ) 

These elements form the matrix of SMFD F a = {/j"}(tv+ 2 )(w+ 2 )- Thus, matrices of 
all layers of T a , L a and F a can be assembled into the corresponding vectors, i.e., Vt = 
{To, Ti, ■ ■ ■ , Vl = { L 0 , Li, • • • , L m _i} and Vp = {F 0 , F x , ■ ■ ■ , F M _ X }. 

Based on the matrices of MTFD, MFPFD and SMFD, many interesting knowledge can 
be found in the open flow systems, such as “trophic levels” and centralities of nodes. 
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TABLE I: A fragment of the NBER-LInited Nations trade dataset 


Year ICode Importer ECode Exporter SITC4 Unit DOT Value Quantity 


2000 457640 Thailand 532760 Germany 

6832 

W 

1 

157 

9 

2000 457640 Thailand 532760 Germany 

6842 

W 

1 

11524 

2731 

2000 532080 Denmark 211240 Canada 

8211 

w 

1 

163 

16 


III. AN EMPIRICAL STUDY: THE MULTI-LAYER FLOW NETWORK OF IN¬ 
TERNATIONAL TRADE 

A. Dataset 

We use the NBER-LInited Nations trade data (http://cid.econ.ucdavis.edu/nberus.html) 
to explore new features of multi-layer flow network of international trade. The dataset covers 
the details of world trade flow from 1962 to 2000, and SITC4 (4-digit Standard International 
Trade Classification, Revision 4) standard [2lJ is used to organize hundreds types of products 
in the dataset. A fragment of the dataset is shown in Table |T[ where “ICode” and “ECode” 
are the corresponding codes for importers and exporters respectively, and “Value” means the 
value of commodity whose unit is thousands of US dollars. “DOT" (direction of trade) has 
two options: 1 (data from the importer) and 2 (data from that exporter). The “Quantity” 
of commodity is measured by “Unit” whose codes can be “W” (weight of metric tons), “V” 
(Volume of cubic meters) and so on. Readers can refer to [19j for a detailed description of 
this dataset. 

In this study, we use the value of commodity as the volume flux of the trade. 

B. Establishment of the multi-layer flow network of international trade 

We use the data in year 2000 to build a multi-layer flow network of international trade. 
In this multi-layer flow network, there are totally 192 nodes (containing two special nodes 
“source” and “sink”) and 1288 different layers, where each common node (and its counter¬ 
parts in the different layers) represents a country and each layer contains a single product’s 
trade flow. Trade flows can be depicted from two different perspectives: the viewpoint of 
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FIG. 2: Illustration of a balanced country node from the viewpoint of commodity flow 


trading commodity flow and that of money flow for trading. A detailed explanation is given 
below. 


1. Viewpoint of trading commodity flow 


In each layer, for the given product of this layer, if there exists a trade relationship 
between two countries, a directed edge can be built from the exporter to the importer, and 
the corresponding value of the trade is set as the volume flux of each edge. The edge from 
“source” to the country node (i.e., j) means the production of this commodity in the country 
j, and the edge from node j to “sink” represents the consumption of this commodity in the 
country j. The volume flux of theses edges are the corresponding value of the production 
(or the consumption). 

For each country node j (except “source” and “sink”) in a certain layer a, according to 
the constraint of balanced inflow and outflow (Fig. [2]), we have the following equation. 

N N 

^2 e ij + e 0j = e % + e l,N +1 ( 9 ) 

i =1 fc=l 

where ■ and e “ Ar+1 are the domestic production and consumption of country j for the 
product in layer a respectively. Because e^- (the flow from the source to j ) and e° N+1 (the 
flow from j to the sink) are not available, for simplicity, they are estimated as below. 


'0 j 


= 0; e' 


j,N +1 


— e Q - e a if e a > e 

Z-U=1 Z-*fc=l c jk 11 1 G? — Z-/fc=l L 


= 0; e oj = Ef= 


-j,N+l 


p a 

1 e jk 


Ef= 


'jk 

e a - 

i 


if e" < V JV 
11 L^i= 1 G? ^ Z^fc= 1 


ij — /-~/k =1 ''jki 
N 


0 a 

' jk■ 
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FIG. 3: Illustration of a balanced country node from the viewpoint of money flow 


2. Viewpoint of money flow for trading 

Besides the viewpoint of trading commodity flow, the flow network of international trade 
also can be established from the viewpoint of money flow. Because the trading of commodity 
is always accompanied by money flow, the direction of money flow is just the opposite of 
that of commodity flow, that is from the importer to the exporter. The edge from “source” 
to the country node j represents country j’s trade deficit driven by the consumption of the 
commodity, and the edge from country node j to “sink” means the surplus of the exports 
over the imports. 

For each balanced country node j (except the source and the sink), Equation [9] still holds 
(Fig. |3j), where e“- and are the money inflow from country i to country j and the money 
outflow from country j to country k respectively, eg) and e^ N+l are the deficit and the 
surplus respectively, which also can be calculated using Equation [lOj 


IV. RESULTS 

A. Trading trophic levels 

1. Countries’ trophic levels from the viewpoints of commodity flow and money flow 

“Trophic level” is a term borrowed from ecology, meaning the position a species occupies 
in a food chain. In flow network systems, we use this variable to indicate the position that 
a node occupied in the whole system, and its value is one node’s distance from “source”. 
Since there may be multiple paths from the source to the node, adopting the distance of the 
shortest path from the source may underestimate the trophic level of a node, and using the 
mean first-passage flow distance (MFPFD) from the source is more reasonable |13j . 


9 


•Honduras 



(a) The viewpoint of commodity flow 


(b) The viewpoint of money flow 


FIG. 4: Trading trophic levels of countries for the commodity of live bovine animals from 

different viewpoints. 


We use the layer of the commodity of live bovine animals as an example. Based on 
MFPFD discussed in the subsection |II B[ we use the MFPFD from the source to each country 
node as its trading trophic level. The trophic levels of each country from the viewpoints of 
commodity flow and money flow are given in Fig J4] respectively. In the figure, country nodes 
are plotted in a circle, whose distances to the center of the circle are proportional to their 
trophic levels, and whose angles and colors are selected randomly. The nodes’ sizes are in 
proportion to the natural logarithm of the nodes’ total outflow (or inflow). 

From the viewpoint of commodity flow, the source of commodity flow is the production 
of this commodity and the sink means the consumption. So for a certain commodity, the 
trophic level of a country node represents the position that the country occupied in the 
global supply chain (i.e., the commodity flow). The smaller the trophic level, the role of 
the country is more inclined to be the producer (i.e., the exporter) of this commodity. On 
the contrary, the bigger the trophic level, the greater the distance from the source to the 
country node and the country is more inclined to be the consumer (that is the importer) 
of this commodity. From Fig|4a| we find that some countries (such as Germany, China, 
Australia, Canada and South Africa), whose trading trophic levels are slightly greater than 
or equal to 1, are inclined to be the exporters of live bovine animals. Some other nodes, such 
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(a) The viewpoint of commodity flow 


(b) The viewpoint of money flow 


FIG. 5: The distribution of trading trophic levels of countries for the commodity of live 

bovine animals 


as Korea Rep., Venezuela and Colombia, which are far away from the centre of the circle, 
are the importers of live bovine animals. The corresponding distribution of trophic levels of 
countries is given in Fig|5a| Obviously, there are three major groups of countries, i.e., the 
countries whose trophic levels are between 1 and 1.6, between 2 and 2.8, and larger than 3. 
The universality of this pattern has been confirmed for a large variety of commodities. 

From the viewpoint of money flow, the source is consumer demand for the commodity 
which leads to trade deficit, and the sink means the surplus. Because money flow can be 
obtained through interchanging “source” and “sink” and reversing the direction of edges of 
commodity flow network, it is intriguing to find some regularities and rules between money 
flow network (say MG a ) and commodity flow network (say CG a ) for the same commodity. 
These regularities and rules include the followings. 

• MG a and CG a have the same MFPFD from the source to the sink (i.e., io,JV+i)- For 
the commodity of live bovine animals, Zq at+i °f MG a and CG a are both 3.3674. 

• For a country node, the larger the trophic level on MG a , the smaller on CG a , and 
vice versa. That is to say, given any two countries i and j, if trophic level of i is larger 
than that of j on MG a , trophic level of i is no bigger than that of j on CG a , and vice 
versa. For instance, Venezuela has a biggest trophic level (which is 3.7134) from the 
viewpoint of commodity flow; on the contrary, its tropic level (which is 1.00) is the 
smallest from the viewpoint of money flow. 
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• The distributions of trade trophic levels of countries from two viewpoints exhibit sim¬ 
ilar patterns: a large percentage of trade trophic levels are equal to or slightly larger 
than 1, and another group of trophic levels is between 2 and 3. 


2. Countries ’ trophic levels for different products 

We compare countries’ trophic levels in different layers with different products. We 
portray them using the countries-products matrix as shown in Figj6} In the figure, each 
point represents a trophic level with the corresponding country and product, and its value is 
depicted by the color, where white indicates the value does not exist, cyan means the value 
is between 1 and 2 (no including 2) and indicating the country tends to be an importer, and 
magenta represents a relatively high value indicating an exporter. Three variables, i.e., the 
number of cyan points (say n cyan ), that of magenta points (say n magenta ) and the mean of 
MFPFD from the source to the country represented by the row (say l qJ, are extracted to 
characterize each row (or column). Rows and columns of the matrix displayed in the figure 
are sorted by the sum of n cya n and n mage nta in descending order, then by 1^. 

In Figj6j we obtain an approximately right angle trapezoidal shape for cyan and magenta 
points. We can found that magenta points are concentrated in the upper portion of the 
trapezoidal shape, while the lower part is mainly cyan dots. It can be interpreted that the 
exporters of a wide spectrum of products are active countries in the international trade, 
which are located in the upper part of the figure and can be regarded as competitive and 
success countries according to [25]. In contrast, bottom rows are those not active in the 
international trading, which can only export a limited variety of products and tend to be 
underdeveloped countries [25]. Economies, which export few types of goods but import a 
large variety of products, also can occur in the upper part of the figure. Examples include 
Saudi Arabia and United Arab Emirates. Besides, we also find some abnormal countries, 
which appear as cyan horizontal lines in the top-right of the figure. It can be interpreted 
that these countries import many kinds of goods, which the vast majority of countries do 
not need to import. It may indicate that these economies are at high risk. A representative 
example is Yugoslavia in the year 2000. 
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FIG. 6: Countries-products matrix reporting countries’ trophic levels in different layers 
with different products. Countries and Products are arranged in descending order of the 
sum of n cyan and n magent a, and then where n cyan , n magen ta and are the number of 
cyan points, that of magenta points and the mean of MFPFD from the source to the 


country represented by the row respectively. 


B. Node centrality 

1. Countries’ node centralities in a certain layer of product 

We calculate country centrality based on flow distances in a certain layer of product. The 
centrality of a country node can be computed as the average of distances from the node to 
all the other countries as given in [13]. However, if one of those distances is infinite, the 
node centrality will become infinite. Thus, the above computation method is infeasible for 
a sparse distance matrix in which most of the elements are infinite. Therefore, we propose 
a new definition of node centrality, called harmonic centrality, as 


N - 1 


( 11 ) 



where ff j (j = 1, ■ ■ ■ , N and j ^ i ) is the SMFDs between nodes i and j, and N is the 
number of common nodes in the flow network. Two special nodes, i.e., source node 0 and 


sink node IV + 1, are excluded, because here we are concerned about the distances of the 


node to the other common nodes, and the distances to the source and the sink have been 


13 













TABLE II: List of countries sorted by fi 


(a) Top 5 countries (b) Bottom 5 countries 


Rank 

Country 

fi 

Rank 

Country 

fi 

1 

Germany 

2.53 

82 

Ecuador 

86.00 

2 

Netherlands 2.74 

83 

Mozambique 86.00 

3 

Hungary 

3.05 

84 

Bahrain 

86.00 

4 

Australia 

3.10 

85 

Qatar 

86.00 

5 

Italy 

3.15 

86 

Singapore 

86.00 


depicted in the metric of trophic level. 

Thus, the harmonic node centrality metric can avoid the infinite distance problem and 
well measure centralities of nodes. A smaller /* implies a more central position the node 
occupied in the flow network, because it has a smaller average distance to the other nodes. 
In Table |TT[ we list the top 5 and the bottom 5 countries with their fi in the decreasing order 
of fi in the money flow network of trading live bovine animals. The top 5 countries are all 
important hub nodes from the viewpoint of topology. 


2. Countries ’ node centralities for different products 

We compute and compare node centralities of countries in different layers of products. 
We establish a countries-products matrix containing the corresponding harmonic centrali¬ 
ties. For example, the row for USA is [3.95,15.50, 5.03, • • • , 7.50], and the column for the 
commodity of live bovine animals is [86.00, n.a., n.a., 7.13, n.a.], where n.a. means 

that the corresponding country do not appear in the trading of the commodity. We then 
transform the matrix into the countries-products matrix containing the rankings of countries 
for each type of commodity. Thus, the row of USA becomes [15, 27,11, • • • ,4], which means 
that USA is ranked 15th, 27th, • • •, and 4th respectively in each column, and the column 
for the product of live bovine animals turns into [79, n.a., n.a., 40, • • •, n.a.]. 

We illustrate the ranking matrix in FigjT], where rows are sorted by the mean of rankings, 
and columns are rearranged in descending order of the number of countries participating 
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Products 

FIG. 7: Countries-products matrix reporting countries’ harmonic centrality rankings in 
columns. Countries are sorted in ascending order by the average ranking, and products are 
arranged in descending order of the number of trading countries. 


in the ranking. From the figure, we get an approximately right angle trapezoidal shape for 
colored dots, which is quite similar to the shape in Figj6j It can be interpreted that active 
countries which trade a wide spectrum of products are likely to have high average harmonic 
centrality; on the contrary, those having few types of trading products tend to have low 
average rankings. Most interestingly, in the figure there are several horizontal color bands, 
which are blue, cyan, yellow and red successively from top to bottom. This phenomenon 
can be named as centrality stratification. It indicates that competitive countries tend to be 
in the center position in the trading of a large variety of products, while underdeveloped 
countries likely rank low in their limited varieties of trading products. We list the top 8 
and bottom 8 economies with their corresponding average harmonic centrality rankings (say 
fi ) and numbers of types of trading products (say n). It implies that fi may be a good 
alternative indicator for countries’ competitiveness from the perspective of node centrality, 
and countries need to promote their harmonic centrality rankings in the trading of various 
products for competitiveness enhancement. 
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TABLE III: List of economies sorted by the average harmonic centrality ranking, n is the 
number of types of trading products and /,- the country’s average harmonic centrality 

ranking. 


(a) Top 8 

Rank 

Economy 

n 

fi 

1 

Germany 

928 

4.95 

2 

USA 

992 

5.47 

3 

France, Monac 946 

6.86 

4 

UK 

938 

7.06 

5 

Italy 

917 

8.04 

6 

China 

914 

11.43 

7 

Spain 

912 

12.13 

8 

Netherlands 

927 

13.01 


(b) Bottom 8 


Rank 

Economy 

n 

fi 

183 

Guinea-Bissau 268 

126.72 

184 

Chad 

129 

126.94 

185 

Greenland 

299 

127.99 

186 

Burundi 

108 

128.60 

187 

Samoa 

123 

133.15 

188 

Occ. Pal. Terr 

61 

134.02 

189 

St.Pierre Mq 

68 

134.88 

190 

CACM NES 

21 

137.86 


C. Mean first-passage flow distances from the source to the sink 


We calculate the mean first-passage flow distances (MFPFD) from the source to the sink 
(say Iq for money flow network in each layer. It indicates the average step that a 
random particle may jump from the source to the sink, and can be regarded as the flow 
length of the network in the layer. In Fig. [8j the distribution and cumulative distribution 
of Iq jv+i f° r a ll 1288 layers are shown. We And that over 27% of Iq N+1 is 3, and the density 
of Iqn+i decreases with the growth of the value of Iq n+1 . The maximal one is 6.24. We 
further compare Iq n+1 among groups of layers containing different categories of products. 
According to the classification of products given by SITC4 [21] , the result is shown in Table 
IV, where the category of 9 (commodities and transactions not classified) is ignored. From 
the table, we And that a group formed by categories 5 (chemicals and related products, 
n.e.s.), 7 (machinery and transport equipment), 6 (manufactured goods classiAed chieAy 
by material) and 8 (miscellaneous manufactured articles) have a signiAcant higher average 
IqN+ i than the group of categories 0 (food and live animals), 1 (beverages and tobacco), 2 
(crude materials, inedible, except fuels), 3 (mineral fuels, lubricants and related materials) 
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FIG. 8: The distribution and cumulative distribution of MFPFDs from the source to the 

sink 


and 4 (animal and vegetable oils, fats and waxes). For simplicity, it can be concluded as 
manufactured products (categories 5, 6, 7 and 8) have a significant larger Iq N+1 than primary 
products (categories 0, 1, 2, 3, 4). This phenomenon also has been confirmed by the stacked 
distribution chart of Z^jv+i f° r different categories of products (Fig. |9j), where blue bars are 
stacked on the left side by contrast with most yellow and red bars located on the right. 


V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we use the approach of flow distances to explore the international trade 
flow system from the perspective of multi-layer flow network. The model of multi-layer flow 
network can well model the open flow system with multiple types of flows, where each layer 
contains one type of flows. We introduce the formal description of multi-layer flow networks 
(MFNs), and flow distances on the MFNs (e.g., the mean first-passage flow distance and 
the mean total flow distance), where a new flow distance called symmetric minimum flow 
distance is proposed. Then, we build the multi-layer flow network of international trade 
from two coupled viewpoints, i.e., the viewpoint of commodity flow and that of money flow. 

Thus, based on the approach of flow distances, some interesting knowledge are discovered 
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TABLE IV: MFPFDs from the source to the sink for different categories of products 


Code Classification 

Average l^ N+1 

Minimal Iq n+1 

Maximal Iq N+1 

0 

Food and live animals 

3.35 

3.00 

5.86 

1 

Beverages and tobacco 

3.31 

3.00 

4.34 

2 

Crude materials, inedible, except fuels 

3.30 

3.00 

4.89 

3 

Mineral fuels, lubricants and related 

3.28 

3.00 

4.80 


materials 




4 

Animal and vegetable oils, fats and 

3.34 

3.00 

4.07 


waxes 




5 

Chemicals and related products, n.e.s. 

3.78 

3.00 

5.61 

6 

Manufactured goods classified chiefly 

3.71 

3.00 

6.17 


by material 




7 

Machinery and transport equipment 

3.75 

3.00 

6.12 

8 

Miscellaneous manufactured articles 

3.63 

3.00 

6.24 

- 

All products 

3.59 

3.00 

6.24 


0.4 

0.35 

0.3 

j? 0.25 

\E 
o 

i 0-2 

Q 0.15 
0.1 

0.05 

0 

2.5 3 3.5 4 4.5 5 5.5 6 6.5 

1 “ 

0,N+1 

FIG. 9: The distribution of MFPFDs from the source to the sink for different categories of 

products 


|0: food and live animals 
1 1: beverages and tobacco 
2: Crude materials, inedible, except fuels 
3: mineral fuels, lubricants and related materials 
4: animal and vegetable oils, fats and waxes 
5: chemicals and related products, n.e.s. 

16: manufactured goods classified chiefly by material 
|7: machinery and transport equipment 
1 8: miscellaneous manufactured articles 
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from the multi-layer flow network of international trade. 

Firstly, countries’ trading “trophic levels” is used to depict the positions that countries 
occupied in the international supply chain. From the distribution of trading “trophic levels”, 
we find that countries can be divided into three groups: countries whose trophic levels are 
slightly bigger than or equal to 1, between 2 and 3, and larger than 3. Besides, since the 
viewpoints of commodity flow and money flow are coupled, some regularities can be found. 
For example, money flow network and commodity flow network for the same commodity 
have the same mean first-passage flow distance from the source to the sink. 

Secondly, by comparing countries’ trophic levels in different layers with different products, 
we find that exporters of a wide spectrum of products are active and competitive countries 
in the international trade, while countries export a limited types of commodities are inactive 
and tend to be underdeveloped countries. Besides, we find some countries import many kinds 
of goods, which the vast majority of countries do not need to import. This phenomenon 
may indicate that these economies are at high economic risk. 

Thirdly, we propose a new node centrality called harmonic centrality for solving the 
problem of infinite distance. A smaller harmonic centrality indicates a more central position 
the node occupied. Then, we compare harmonic centralities of countries in different layers of 
products. It is interesting to find the phenomenon of centrality stratification. It means that 
competitive countries tend to be in the center position in the trading of a large variety of 
products, while underdeveloped countries likely rank low in their limited varieties of trading 
products. 

Fourthly, we compute the mean first-passage flow distances from the source to the sink 
for different types of commodities, which can be regarded as the flow length of the network 
in each layer. We find that manufactured products have significant larger mean first-passage 
flow distances from the source to the sink than primary products. 

Our findings demonstrate the effectiveness of the proposed model of multi-layer flow 
networks and the approach of flow distances. 
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